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We study the dynamics of an exactly solvable lattice model for inhomogeneous 
interface growth. The interface grows deterministically with constant velocity 
except along a defect line where the growth process is random. We obtain exact 
expressions for the average height and height fluctuations as functions of space and 
time for an initially flat interface. For a given defect strength there is a critical 
angle between the defect line and the growth direction above which a cusp in 
the interface develops. In the mapping to polymers in random media this is an 
example for the transverse Meissner effect. Fluctuations around the mean shape 
of the interface are Gaussian. 
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The problem of nonlinear, KPZ-type interface growth in 1+1 dimensions 
and the closely related problem of polymers in random media has been the 
subject of many studies over the past decade |]]. The underlying assumption 
in this approach is that growth occurs stochastically normal to the inter- 
face. The dynamics of the interface are then given by the KPZ equation 
dth = uV 2 h + X(Vh) 2 + 77 0. The deterministic version of the equation 
(without the noise 77) can be solved explicitly Among other things one finds 
that in the infinite time limit the interface becomes flat. However, in the 
presence of noise, even an initially flat interface roughens and the deriva- 
tion of quantities related to the roughening process and other fluctuation 
phenomena become an interesting and challenging problem, usually tackled 
using either the renormalization group or the study of lattice models. In this 
context the asymmetric exclusion process || , a driven lattice gas model, has 
played a special role. Both numerical and exact analytical results have been 
obtained from this microscopically motivated exactly solvable model, which 
maps to a lattice growth model in the universality class of the noisy KPZ 
equation. 

Recently Kallabis and Lassig studied the KPZ equation with spatially 
localized noise and found an interesting phase diagram even in one dimen- 
sion [|J. While this work focussed on the properties of the directed polymer 
mapping, subsequent work elucidated some properties of the growth model 
itself ||. Among other things it was found that in one dimension the inter- 
face in the steady state has a cusp with constant slope for arbitrarily small 
noise. This observation is, in fact, not new: Already some while ago, an 
exclusion process on a lattice was introduced which mimics the situation of 
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homogeneous and deterministic growth except along one defect line where 
growth is stochastic 0. In the steady state of this non-equilibrium system 
one finds indeed such a cusp, again for arbitrarily small noise if the defect 
line is parallel to the average growth direction as assumed in the continuum 
model studied in Refs. [|J [5] . This makes a study of the dynamical properties 
of the lattice model desirable. It is the aim of this paper to study the time 
evolution of an initially flat interface using this model. 

The model introduced in Ref. || is an exactly solvable version of the 
deterministic asymmetric exclusion process on a ring with a defect line. We 
use the well-known RSOS mapping from an interface to a lattice gas |7], §]: 
In the interface model the (integer) height variables h(y, t) may differ on 
neighbouring (integer) sites y, y + 1 only by ±1. As a result of this restriction, 
growth can occur only in local minima, no overhangs can develop. If a growth 
event occurs, the height at site y increases by two units (see Fig. 1). In lattice 
gas language the height differences between neighbouring sites are mapped to 
a particle occupation number n x (t) = 0, 1 with the presence of a particle on 
x corresponding to slope -1 between sites y — 1 and y in the interface model 
and a vacancy at site x corresponding to slope +1 (see Fig. 1). Growth at 
site y in the growth model then corresponds to a particle hopping from site 
x to a vacant site x + 1. Note that each lattice site can be occupied by at 
most one particle. In a finite system of L = 2N sites for the particle model 
one needs N particles in order to ensure periodic boundary conditions in the 
height model. In an infinite system, an average density p ^ 1/2 over some 
finite range V corresponds to a non-zero average slope 1 — 2p in that range. 

So far this mapping keeps track only of the height differences between 
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neighbouring sites. In order to make the two models equivalent, one has 
to introduce an additional integer random variable ho G Z in the particle 
system which gives the absolute height of the interface at some arbitrary, 
but specified point, e.g. y = 0. In the lattice gas mapping the value of this 
random variable is increased by two units each time a particle hops from 
site to site 1. From this one may reproduce the height at any point since 
h(y) = h + E v x =i(l-2n x ). 

The deterministic time evolution is realized by a parallel updating scheme 
in which in a first step all odd pairs of sites (2x — 1, 2x) are updated according 
to the following rules: If there is a particle on site 2x — l and a vacancy on site 
2x, then the particles hop with probability 1 to site 2x. If the pair of sites is in 
one of the three remaining configurations, nothing changes. These rules are 
applied in parallel to all such pairs. In a next step one shifts the pairing by 
one lattice unit and applies the same rules to the even pairs (2x, 2x+l). This 
completes one full time step. Note that so far this model is fully determin- 
istic, and any initial configuration will eventually evolve into a flat interface 
characterized by a particle configuration (. . . , 1, 0, 1, 0, 1, 0, 1,0,...). The dy- 
namical properties of a similar system with slightly different deterministic 
updating rules were studied in ||. 

In order to introduce a defect line at y = we change the rules such 
that a particle on site x = hops only with probability p < 1, if site 1 
is vacant. This randomness results in non-trivial behaviour of the system. 
The stationary properties of this model on a finite ring with N particles 
were studied in || (JV arbitrary). Here we shall consider the dynamics of 
the infinite system with flat initial configuration (...,1,0, 1, 0, 1, 0, 1,0,.. .) 
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where the particles are on the odd sites while the vacancies are on the even 
sites. Because of reflection symmetry one has h(y) = h(—y). 

It is useful to note that the time evolution may be written in terms of 
a transfer matrix T = pTo + (1 — p)T\ acting on states corresponding to 
some given configuration of the system ||. Any state is fully character- 
ized by the particle occupation numbers n x and the integer height variable 
ho- Tq is the transfer matrix for the system without defect, while T\ is the 
transfer matrix of the system with full blockage where a particle at site 
can never move to site 1, even if it is vacant. Acting with To on the state 
(. . . , 1, 0, 1, 0, 1, 0, 1,0,.. .) reproduces this state by shifting the whole con- 
figuration by two lattice units to the right. This is in accordance with the 
remark that this configuration is stationary if p = 1. The height variable ho 
increases by two units since a particle went from site two site 1. This means 
that the interface has grown everywhere by two units. On the other hand, the 
action of T\ on this state results in a configuration (. . . , 0, 0, 1, 0, 1,0,1,0.. .) 
where the first occupation number in this set (which is 0) corresponds to the 
occupation on site 1. The occupation numbers on site and on the negative 
sites follow from symmetry. In this case ho has not increased, only h(y) for 
y < 1 has grown by two units. Taking the t th power of T yields all possible 
configurations the interface may take after t time steps. Each realization 
acquires a factor p k (l —p) l ~ k which gives the probability that this particular 
growth history has taken place. 

Using the decomposition of T into T and 7\ it becomes straightforward 
to calculate to respective probabilities for all these configurations: (1) The 
odd sublattice remains empty for all times. (2) The consequence of the 
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defect for the even positive sublattice is the injection of a particle on site 
with probability p in each time step. (3) If a particle has been injected 
it will move in each subsequent time step by two lattice units without any 
interaction with other particles. As an example for what this corresponds to 
in height language, consider ho which we take to be zero at time t = 0. With 
the dynamical rules as explained above one finds after one step ho = 1 with 
probability p and ho = with probability 1 — p. After two steps one finds 
ho = 2 with probability p 2 , ho = 1 with probability 2p(l — p) and ho = 
with probability (1 — p) 2 , and so on. By construction, the height on the even 
sublattice will always be even, while on the odd sublattice it will always be 
odd. Let us denote h(2y,t) = 2h{2y,t). Then one gets for the probability 
P y ,t{h) of finding the height h at site y at time t (for y > 0) 



On the odd (positive) lattice sites one has always, i.e. for all times and all 
realizations of the randomness, h(2y — l,t) = h(2y,t) — 1. For the negative 
half space y < one gets P[h(—y,t)] = P[h(y,t)] by symmetry. This result 
gives a complete description of the time evolution of the local heights of an 
interface which is initially (macroscopically) flat.| Note that there are no 
height fluctuations at all for y > t. This must be so since the perturbation 
of the interface caused by the defect spreads with finite velocity v = 1 which 
is the bulk growth velocity. 

2 Eq. (1) is easy to verify e.g. by applying T two or three times to the initial state. 
This makes the structure of the interface after t steps already quite clear. 




k,t (v>t) 



(i) 
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In what follows we study only the non-negative even sublattice up to site 
y = t, since all quantities relating the odd sublattice, the negative half space 
and the region y > t are then trivially given. From the distribution ([I]) one 
obtains the generating function for height fluctuations 

(e ah ^) = e 2ay {l-p + pe 2a y- y (0 < y < t) 

= e 2at (y>t) (2) 

and in particular the average height 

(h(2y,t)} = 2t - 2(1 -p)(t-y) (0 < y < t) 

= 2t (y>t). (3) 

The height fluctuations can be obtained by taking the second logarithmic 
derivative of the generating function. One finds 

(h 2 (2y,t))-(h(2y,t)) 2 = Ap(l-p)(t-y) (0 < y < t) 

= (y>t) (4) 

From Eq. (|3|) one finds that the defect causes the development of a cusp with 
constant slope 2(1— p). At site y = the interface grows with constant speed 
Vq = 2p, while in the bulk, at distances larger than t, it grows with velocity 
v b = 2 (Fig. 2). From ([]]) one finds that the fluctuations of the interface 
round its mean value @ are Gaussian (for t large) with a space-dependent 
variance given by (|j). These are the main results of this paper. 

We conclude by briefly discussing the situation in which the defect line 
is not parallel to the growth direction of the interface. This can be realized 
within this model by taking a (macroscopically) non-zero initial slope s = 
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1 — 2p of the interface. Since growth occurs normal to the interface the net 
result is indeed a defect line tilted relative to the main growth direction. In 
particle language one has to choose an initial state with a density p 7^ 1/2., 
e.g. a configuration (...,1,0, 0, 0, 1, 0, 0, 0, . . .). From the analysis of the 
steady state of this system one finds that a phase transition takes place at 
p c = 2p for p < 1/2 (positive initial slope) and p c = 2(1 — p) for p > 1/2 
(negative initial slope) respectively || . For the dynamics of an initially flat 
interface the steady state properties of the system suggest that for p > p c 
the interface remains flat for all times except in a finite region close to the 
defect. For p < p c a cusp will evolve. Unlike in the situation described above, 
the boundary of this cusp will not be sharp, but fluctuate. In a finite system 
these fluctuations scale in the system size like L 1 / 2 and one would therefore 
expect a power growth t 1//2 in time of these domain boundary fluctuations in 
the infinite system. For fixed p this is a phase transition that takes place at 
critical slopes s c = ±(1 — p). For — (1 — p) < s < (1 — p) the system develops 
the cusp. 

This situation is also of interest in the well-known mapping of this prob- 
lem to a directed polymer in a medium with a random defect line, as studied 
in Ref. || for the special case p = 1/2 where the defect line and the av- 
erage direction of the polymer are parallel. The development of the cusp 
corresponds to having a bound state where the polymer is bound to the (at- 
tractive) defect line. It is no surprise that this happens at arbitrarily small 
defect strength. This simply reflects the fact that in one dimension a delta 
function potential (the defect line) with arbitrarily small amplitude has a 
bound state. Translated into polymer language the phase transition in the 
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interface model is a transition from an bound state of the polymer for a tilt 
of the defect line against the (average) direction of the polymer below the 
critical angle s c to an unbound state for a tilt above s c . This is a simple 
example for the occurrence of the transverse Meissner effect |TU . 



A macroscopic cusp with constant slope appears to be an universal feature 
of growth models with a defect line where the growth process is impeded. In 
this sense some features of the model studied in this paper may be of more 
than purely theoretical interest. Normal growth, which is the underlying 
assumption that leads to the KPZ equation, has been reported to occur e.g. 
in experiments on wetting of tapes plunged in a liquid JO] . It would therefore 
be interesting to adapt the model to this situation by an appropriate choice 
of boundary conditions. On the other hand, our simple model has some clear 
limitations. The strict discontinuity of the slope of the interface at y = t (for 
the half-filled system) is obviously an artefact of the dynamics chosen here 
and cannot be expected to be found in any real system. One would rather 
expect a smeared out discontinuity as in the system with a tilted defect or 
as in the fully stochastic version of the model with a defect |H| . Another, 
intrinsic, limitation of the model is that growth at the defect is always slower 
than in the bulk. It would be very interesting to study a model where the 
defect growth could be made faster than the bulk growth velocity. 

The author would like to thank M. Lassig and H. Kallabis for useful 
discussions. This work was supported by an EC Fellowship under the Human 
Capital and Mobility program. 
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List of Figure Captions 



Fig. 1: The mapping between the restricted interface and the particle exclu- 
sion process: we show a possible interface configuration and the correspond- 
ing particle occupancies on a lattice with sites labeled by y. The indicated 
flips in the interface correspond to particles hopping on the lattice, marked 
by horizontal arrows. 



Fig. 2: Macroscopic average interface shape at time t — (flat interface with 
height 0) and at some later time t. The position of the defect at y = is 
indicated by the vertical dashed line. The defect line is parallel to the bulk 
growth direction. 
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Figure 2: 
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